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Abstract. We present a quantum Monte Carlo study of two-dimensional dipolar Bose gases in the limit of zero temper- 
ature. The analysis is mainly focused on the anisotropy effects induced in the homogeneous gas when the polarization 
angle with respect to the plane is changed. We restrict our study to the regime where the dipolar interaction is strictly 
repulsive, although the strength of the pair repulsion depends on the vector interparticle distance. Our results show that 
the effect of the anisotropy in the energy per particle scales with the gas parameter at low densities as expected, and 
that this scaling is preserved for all polarization angles even at the largest densities considered here. We also evaluate 
the excitation spectrum of the dipolar Bose gas in the context of the Feynman approximation and compare the results 
obtained with the Bogoliubov ones. As expected, we find that these two approximations agree at very low densities, 
while they start to deviate from each other as the density increases. For the largest densities studied, we observe a 
significant influence of the anisotropy of the dipole-dipole interaction in the excitation spectrum. 
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1 Introduction 

In 2005 Griesmaier et al JJ] and Stuhler et al JUl obtained the 
first Bose-Einstein condensation of 52 Cr atoms, an achievement 
that triggered the advent of a lot of new theoretical and exper- 
imental work on dipolar Bose-Einstein condensates (DBEC). 
The main difference between 52 Cr condensates and previous 
ones with other alkalis is the moderately large magnetic dipolar 
momentum of Chromium (6/Ib) which makes dipolar interac- 
tion comparable in strength to the usual Van der Waals forces, 
thus producing significant new effects that can be measured in 
experiments. 

Dipolar quantum gases are interesting from the experimen- 
tal and theoretical points of view due to the novel features that 
the dipolar interaction introduces: anisotropy and long range. 
Dipolar interaction can be either attractive or repulsive depend- 
ing on the relative orientation of both the position and the dipo- 
lar moments of the particles, a fact that can crucially affect the 
stability of the many-body system and introduces a new degree 
of freedom which enriches the phase diagram. On the other 
hand, the long range character of the dipolar interaction leads 
to scattering properties that are radically different from those 
found on the usual short-ranged potentials of quantum gases. 

On the experimental side there are new exciting results com- 
ing from the work with polar molecules on one sideQHI!] 
and on Dysprosium la, 01 and Erbium [8] condensates on the 
other. Polar molecules have a large and tunable electric dipole 
moment, but are difficult to cool down to quantum degeneracy 
while keeping the system stable. On the other hand, Dyspro- 
sium and Erbium have magnetic moments comparable to 52 Cr 



but a significantly larger mass, leading to a much larger dipolar 
coupling constant with the added benefits of producing a neat 
and stable magnetic dipolar condensate where dipolar effects 
are more important than the Van der Waals forces. 

Dipolar quantum gases are also challenging from the theo- 
retical point of view. Many works studying several character- 
istics of these systems have been published in the last years. It 
is known that an homogeneous dipolar quantum gas is dynam- 
ically unstable against collapse in three dimensions, while the 
trapped case is conditionally stable depending on the geometry 
of the trapping potential l9l [loll . This implies that the DBEC 
stability is enhanced in pancake-like traps, where the dipolar 
interaction is globally repulsive. Low-dimensional dipolar gases 
have also gathered major theoretical interest because the typi- 
cal experimental setup involves strongly anisotropic traps to 
stabilize the system. In particular, these trapping potentials can 
be tight enough to make the system effectively two- or one- 
dimensional. Many interesting studies concerning two- or quasi- 
two dimensional dipolar quantum gases have been performed 
in recent y ears, inclu ding the analysis of two-body scattering 
properties 
system H 



HUTU Hit], static properties of the many body trapped 
and homogeneous gas lfl5l[Tr^[T^[l^[l9ll . and 
some works about the dynamic response Il20l I2U 12211 . 

In this work, we analyze the two-dimensional quantum gas 
of fully polarized bosonic dipoles. The dipolar interaction, V c / C / (r), 
between two dipoles is given by 



Pi-p2-3(pi-f)(p 2 -f) 



(1) 
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with r the relative position vector between them and Cdd the 
constant defining the strength of the dipolar interaction. For 
permanent magnetic dipoles, Cdd = £ioA* 2 where Hq is the per- 
meability of vacuum and fi is the dipole moment of the atoms. 
Alternatively, the electric dipole moment can be induced by 
an electric field E, and in this case the coupling constant is 
Cdd — d 2 /fit), where d = 6cE with a the static polarizability 
and £o the permittivity of vacuum. For a system of fully po- 
larized dipoles in 2D as the ones considered in this work, P[ 
and p 2 are parallel to the polarization field (which lays on the 
xz -plane), and form an angle a with the normal direction to the 
plane, defining a fixed direction in space, see Fig.Q] In this case 
Vddi r ) simplifies to 



C dd l-3A 2 cos 2 
An r 3 



(2) 



where A = sin a and (r, 6) the in-plane distance and polar an- 
gle, respectively. Notice that, since a is fixed, A < 1 is a con- 
stant of the problem. As usual in the study of dipolar gases, we 
express all quantities in dipolar units, dipolar length and en- 
ergy, which are given by r c / = mCdd /ATth 2 and Ed = h 2 /mr 2 l , 
respectively. Another important feature of the dipolar interac- 
tion in two dimensions is that, contrarily to the 3D case, it is 
short ranged. We restrict our study to polarization angles where 
the interaction is fully repulsive, i. e. a < a c w 0.615, where 
the stability of the system against collapse is ensured but the 
interaction may still show a high degree of anisotropy. 

By means of diffusion Monte Carlo (DMC) simulations we 
evaluate the equation of state of the system, extending previous 
results lfl9ll up to values of the density well above the mean- 
field regime, and show the effects of the anisotropy in the en- 
ergy per particle and in some of the most relevant ground-state 
structural quantities of the system. We also present the excita- 
tion spectrum of the anisotropic gas in both the Bogoliubov and 
Feynman [23] approximations. We compare these two approx- 
imations at low density, showing that they both agree well with 
each other. At higher densities, the Bogoliubov approximation 
breaks down, so we use the Feynman approximation to show 
how a roton minimum emerges and develops differently as a 
function of the direction in momentum space. 

The rest of the paper is organized as follows. In Section II, 
we present DMC results for the ground state of the 2D gas of 
tilted dipoles. The effects of the anisotropic interaction in both 
the energy and structure properties when the density increases 
are discussed. In Section III, we calculate the excitation spec- 
trum using the Feynman approximation, relying on the DMC 
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Fig. 1. (Color online) Left: Two dipoles confined to move on the X-Y 
plane. The polarization field lays on the XZ plane and fixes a direction 
in space forming an angle a with the z axis, r and 9 are polar coor- 
dinates in the X-Y plane. Right: Particular case of two dipoles on the 
X-Z plane (y = 0). 



results for the static structure functions, and compare it with the 
Bogoliubov spectrum. Also in this case, relevant signatures of 
the anisotropy are observed, mainly around the roton momenta. 
Finally, in Section IV we present the summary and conclusions 
of our work. 



2 Ground state: energy and structure 

We have studied the many -body properties of a two-dimensional 
gas of bosonic dipoles using the DMC method. DMC is a zero- 
temperature first-principles stochastic method which leads to 
exact properties of the ground state of bosonic systems. It is a 
form of Green's Function Monte Carlo which samples the pro- 
jection of the ground state from the initial configuration with 
the operator exp[— (J*f — £o)t]. Here, ^ is the Hamiltonian 
of the system, Eq is a norm-preserving adjustable constant, and 
T is the variable which corresponds to imaginary time. The 
simulation is performed by evolving in time T by means of a 
combination of diffusion, drift and branching steps acting on 
walkers (sets of 2N coordinates) representing the wavefunction 
of the system. 

The Hamiltonian of of N fully polarized dipoles in two di- 
mensions, written in dipolar units, is given by 



1 N 



I 



1 -3 A 2 cos 2 ft. 



(3) 



where ry and ft ; are the relative distance and polar angle formed 
by dipoles i and j, respectively. This Hamiltonian is valid only 
when Van der Waals interaction can be neglected in front of 
dipolar forces. As we have commented in the previous Section, 
we restrict our study to polarization angles a < a c = 0.615, 
thus ensuring the potential is fully repulsive and the system is 
stable. 

In order to guide properly the diffusion process and im- 
prove the variance of the results one introduces in the DMC 
method a trial wave function for importance sampling. In the 
present case, we use a variational Jastrow wave function of the 
form 

y(ri,...,r*)=n/(ry) ) (4) 

where ry = r, — r,. Notice that due to the anisotropy of the 
system, the two-body correlation factor /(r) depends on the 
whole vector, rather than on its magnitude only. If the den- 
sity is not high, the low-energy two-body scattering solution 
greatly influences the properties of the many body system, as 
three body scattering processes have extremely low probabil- 
ity. For this reason, we use as a Jastrow two-body correlation 
function the anisotropic zero-energy scattering solution fl9ll 
matched at some large distance, £, , with the symmetrized form 
of a phononic wave function fph(/ij) = exp (—CI ry) . l24ll B, be- 
ing a variational parameter. 

Figure|2]shows the ratio of the energy per particle obtained 
from our DMC calculations and the mean field prediction ll25ll 
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Fig. 2. (Color online) Ratio of the energy per particle of the gas of 
dipoles for different polarization angles to the mean field prediction 
((5) as a function of the gas parameter x = na 2 . Red circles, blue trian- 
gles and black inverted triangles correspond to a = 0.2,0.4 and 0.6, 
respectively. The line corresponds to the mean field prediction. 
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Fig. 3. (Color online) Static structure function for polarization angles 



a = 0.4 and a 



■ 0.6 for nr 2 . 



2,2 and 2 . Red squares and blue 



circles show the two cuts S(k,0) and 5(0, k), respectively. 



where x = na 2 is the gas parameter, with n and a being the den- 
sity and the s-wave scattering length, respectively. We report 
results for three polarization angles a = 0.2,0.4,0.6. It is im- 
portant to notice that, for a given value of the gas parameter x, 
different polarization angles imply different scattering length 
values and different densities. However, and as it can be seen 
from the figure, all energies corresponding to the same x col- 
lapse into the same curve with very small deviations even at 
the higher densities considered. That means that the effects of 
the anisotropy of the interaction are accurately contained in the 
polarization dependent scattering length a(X), which is well 
approximated by the law lfl9ll 



(6) 



where y = 0.5772... is the Euler-Mascheroni constant and the 
factor e 2 ? corresponds to the s-wave scattering length of the 
isotropic dipolar system in 2D (a = 0) ifTTll . The scaling law 
observed in the energy, as reported in Fig. [2] may be also at- 
tributed to the small contributions to the energy coming from 
anisotropic terms (explicit 9 terms) contained in the local en- 
ergy estimator. 

The anisotropic character of the dipole-dipole interaction 
has a direct influence on the ground-state wave function that 
is mirrored in the ground state expectation values of many- 
body operators. DMC allows us to evaluate pure estimations 
l26ll of these observables. Contrarily to isotropic fluids, the 
static structure factor S(k), in the dipolar gas depends on the 
full vector k = (k x ,k y ) rather than on its magnitude. Figure [3] 
shows two cuts of 5(k) along the perpendicular and parallel 
directions with respect to the polarization plane, correspond- 
ing to the lines where the interaction is most and least repul- 
sive, respectively. As expected, the effect of the anisotropy is 
more clearly seen at high densities and large polarization an- 
gles. In particular, for fixed a the separation between S(k,Q) 
and 5(0, k) is enhanced with increasing density. In much the 
same way the separation between the two cuts of the structure 
factor also increases for fixed density when the polarization an- 
gle is increased. At the largest densities considered, the effect 
of increasing the polarization angle makes the peak in 5(k) in- 
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Fig. 4. (Color online) Two-body distribution function for polarization 
angles a = 0.4 and a = 0.6 for nr 2 , = 2 4 ,2 and 2°. Red squares and 
blue circles show the two cuts g(r,0) and g(0, r), respectively. 



crease in the strongly interacting direction and decrease in the 
weakly interacting one. This behavior is clearly observed in the 
Figure for the case nr 2 d = 2 4 and a = 0.6. 

The dependence of the interaction potential on the polar 
angle 9 induces also changes in the distribution functions in 
different directions. We have calculated the two-body distri- 
bution function g(r) of the gas for different densities and po- 
larization angles to visualize the effects of anisotropy in the 
spatial structure. As in the case of the static structure function, 
we have selected two cuts, g(r,0) and g(0,r), corresponding to 
the less and more repulsive directions. Selected results for the 
same densities and angles reported in Fig. [3] are shown in Fig. 
|4] As one can see, the anisotropy is also observed in the spa- 
tial structure but the effect is smaller than the one observed in 
S(k). Nevertheless, looking at the case (nr 2 j = 2 4 ,a = 0.6) one 
observes that g(r,0) has less structure than g(Q,r) which has a 
well defined first peak. 
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3 Excitation spectrum 

A relevant issue in the study of tilted dipolar gases is the influ- 
ence of the anisotropy of the interaction on the collective exci- 
tation spectrum. In this Section, we analyze this problem within 
two standard methods used currently in the study of Bose flu- 
ids: the Feynman and Bogoliubov approximations. 

The Feynman spectrum is easy to derive from a simple sum 
rules argument and provides a single line in (k, co) space corre- 
sponding to a set of infinite lifetime excitations of energy l23ll 
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In this approximation, e(k) depends directly on the static struc- 
ture factor, the only non-trivial quantity, a nd p rovides an upper 
bound to the actual excitation spectrum ll27ll . In systems like 
liquid 4 He, this bound is closer to the experimental mode the 
lower the total momentum is. 

On the other hand, we can study the excitation spectrum of 
the low density two-dimensional dipolar gas in the framework 
of the mean-field theory using the 2D time-dependent Gross- 
Pitaevskii equation, 



^dw Ti 2 ~ 



4xtr 



(8) 



25J. Per- 



Fig. 5. (Color online) Comparison of Feynman (symbols) and Bogoli- 
ubov (lines) excitation spectrum for angles a = 0.4 and a = 0.6 and 
nr^ = 2~ 4 ,2~ 6 and 2~ 8 . Red solid and blue dashed curves show the 
two cuts e(k,0) and e(0,k) corresponding to Bogoliubov approxima- 
tion respectively. Red open squares and Blue open circles show e(k, 0) 
and £(0,k) corresponding to Feynman approximation respectively 



The functional form of the pseudopotential V ps (r) ( fTOb as a sum 
of two terms, one isotropic and another anisotropic, follows the 
same prescription used in the three-dimensional problem l28l 
l29tl . One can consider a linear perturbation of the condensate 
wave function of the system of the form 



where g is the 2D coupling constant g 

r"6 — 

forming a standard Bogoliubov-deGennes linearization one finds 
the well-known Bogoiubov spectrum 



e(k) 



lh 2 k 2 fh 2 k 2 



2m V 2m 



■2gn 



(9) 



\ff(x,t) = e ^' (y/n+ Sy(x,t)) , 
where the perturbative term 8y/(x,t) is given by 
5y(V) = ce' (kx - fl,f) , 



Although the spectrum obtained using this approach has con- 
tributions coming from the anisotropic character of the inter- 
action due to the polarization angle dependence of the scatter- 
ing length, not all contributions of the same order are taken 
into account. This simple Bogoliubov approach disregards the 
contribution coming from higher angular momentum channels, 
keeping only s-wave scattering processes. However, we know 
that different angular momentum channels couple in a non- 
trivial way in a dipolar system and so we have to take them 
into account. We know from the analysis of the zero-energy 
two-body problem that higher order partial wave contributions 
appear with higher orders in A 2 , so the leading corrections ap- 
pear in c/-wave. In order to consider the contribution of the d- 
wave we use the following pseudo-potential 



y ps (r)=^ 2 )(r)- 



Cdd 3 A" cos 29 
An r 3 



(10) 



that leads to the following Gross-Pitaevskii equation 

dV' ft 2 «-,2 , ,2 



Cdd 

4tz 



,3A 2 cos20 , 2 



v>\3 



(ID 



(12) 



(13) 



where c is the (small) perturbation amplitude. 

By inserting (fT2l into Eq. (fTTT i. and neglecting non-linear 
terms, one finds the equation fulfilled by the small perturbation 
8\jf, 



ih 



d8\jf 



n 2 o 

= - —V z 8\ir+(2gn- Ii)8\j/ + gn8\jf* 
2m 

n{F(k)8w*+F(-k)8xir) , 



(14) 



where F(k) is given by 

Cdd 



F(k) 



3A 2 cos20 ik 



(15) 



and y = r r'. Now, taking into account that for a dilute homo- 
geneous system the chemical potential is ji = gn and that for 
the two-dimensional dipole-dipole interaction F(k) = F(— k), 
we finally arrive at the following expression for the Bogoliubov 
spectrum 



h 2 k 2 fh 2 k 2 

£ ( k ) = \Hzr\ + ^A 2 cos20 k ) 



2m V 2m 



(16) 



where 9^ is the angle formed by the momentum of the excita- 
tion and the jt-axis. 
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Fig. 6. (Color online) Feynman excitation spectrum for angles a = 0.4 
and a = 0.6 for nr^ = 2 4 ,2 and 2°. The red open squares and blue 
open circles show the two cuts e(fc,0) and e(0,k), respectively. 



The comparison between the Bogoliubov approximation gi- 
ven in this expression and the excitation spectrum obtained 
from DMC calculations using the Feynman approximation is 
shown in Fig. for several values of the density and polar- 
ization angle. We can see from Fig. [5] that, as expected, the 
Bogoliubov and Feynman approximations coincide at very low 
densities. It is also noticeable the fact that, for a given value 
of the density, the Bogoliubov approximation is closer to the 
Feynman prediction at large polarization angles. This stresses 
once again that the relevant quantity describing the low density 
dipolar Bose gas is the gas parameter x = na 2 , and according to 
equation (O the scattering length a(X) decreases with increas- 
ing polarization angle. For a fixed density, x decreases when a 
increases, and the Feynman prediction gets closer to the Bo- 
goliubov mode, which is know to successfully characterize the 
excitation spectrum of Bose gases when x — >• 0. We can con- 
clude from Fig. [5] that Feynman and Bogoliubov approxima- 
tions are close to each other at small values of the momentum 
k. Finally, one also sees that the excitation spectrum becomes 
isotropic when k — > 0. 

Furthermore, the Bogoliubov approximation is expected to 
be valid only at very low densities while the Feynman approx- 
imation is known to provide an upper bound to the exact ex- 
citation spectrum of the system. To have some insight on how 
e(k) evolves with the density we show in Figure[6]the Feynman 
mode at higher values of n. The results presented in the figure 
correspond to densities that are still far from the crystallization 
point of the isotropic system lfl5l [l6ll . From Fig. |6l one can 
see that with increasing density the spectrum develops a roton- 
like minimum which for fixed density and polarization angle is 
deeper in the most repulsive direction. It is interesting to notice 
that as the anisotropy of the interaction is increased,i.e., when 
the polarization angle grows, the roton minimum is deeper in 
the more repulsive direction while in the orthogonal direction 
the spectrum does not show any minimum in the range of con- 
sidered densities. In fact, the emergence of the roton and its 
eventual zero-energy limit has been discussed as a clear signa- 
ture of the instability of the system when the critical polariza 
tion angle is higher than a, 
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4 Summary and Conclusions 

To summarize, in this work we have described some proper- 
ties of a two-dimensional dipolar Bose gas where the polar- 
ization field forms an angle a with the normal direction. The 
projection of the polarization vector on the plane defines the 
x-axis, where the interaction potential is weaker than in any 
other direction. In this situation there is a critical polarization 
angle a c = 0.615 where the potential starts to show attractive 
regions. We have used the two-body zero energy wave function 
to build a many-body Bijl-Jastrow wave function that we used 
as an input for our diffusion Monte Carlo simulations of the 
homogeneous polarized dipolar Bose gas. We have presented 
results for the energy per particle in comparison with the well 
known mean field prediction of the two dimensional equation 
of state at low densities. The scaling of the energy in the gas 
parameter is preserved up to values of x much larger than those 
of the mean-field regime, corresponding to values x < 10 -2 , in 
which the energy only depends on the gas parameter and not 
on the specific interaction. 

In the second part, we have studied the excitation spectrum 
of the system in two different frameworks, the Bogoliubov and 
Feynman approximations. We have shown that the agreement 
between these two approaches is very good at very low den- 
sities as it was expected. We have derived the first anisotro- 
pic correction to the Bogoliubov spectrum, and we have shown 
that in the range of validity of the Bogoliubov approximation 
the anisotropy plays an extremely small role. At higher densi- 
ties, the Feynman approximation predicts a strongly anisotro- 
pic spectrum showing a roton-like minimum that depends on 
the direction of the momentum. 
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